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Abstract —The main usage of cellular networks has changed 
from voice to data traffic, mostly requested hy static users. In 
this paper, we analyze how a cellular network should he designed 
to provide such wireless hroadhand access with maximal energy 
efficiency (EE). Using stochastic geometry and a detailed power 
consumption model, we optimize the density of access points 
(APs), number of antennas and users per AP, and transmission 
power for maximal EE. Small cells are of course a key technology 
in this direction, hut the analysis shows that the EE improvement 
of a small-cell network saturates quickly with the AP density and 
then “massive MIMO” techniques can further improve the EE. 

I. Introduction 

The data traffic in cellular networks has experienced an 
exponential growth in the past couple of decades, and this 
trend is expected to continue in the foreseeable future ||T]. 
Whenever one observes an exponential growth rate, one has 
to question whether it can be sustained—in particular from 
an ecological perspective. The energy consumption of the 
information and communication technology (ICT) industry 
and its related pollution have already become major societal 
and economical concerns 13. To accommodate a 1000 x 
higher data traffic over the next 10-15 years, without increasing 
the ICT footprint, we need to design new technologies that 
improve the overall energy efficiency (EE) by 1000 x Q- 

Important steps towards improving the EE of cellular net¬ 
work were taken in 0-l^> where a methodology for mea¬ 
suring the energy consumption was developed. It shows that 
the radiated signals, radio-frequency (RE) circuits, baseband 
processing, and backhaul infrastructure are all contributing to 
the overall energy consumption. A promising way to improve 
the EE is the small-cell network approach ||7|, which is based 
on deploying a high density of low-power access points (APs) 
endowed with multiple antennas and advanced beamforming 
capabilities. Smaller cells require less radiated signal energy 
since the distances between APs and user equipments (UEs) 
are shorter. This comes at the price of deploying many more 
APs, which increase the circuit energy consumption instead 
Another promising approach to improve the EE is massive 
MIMO (multiple-input multiple-output) 0’ where 

each AP serves tens of UEs in parallel by beamforming from 
hundreds of small antennas. The form factor of the antenna 
array is not “massive” but may be as a flat-screen television or 
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smaller. The array gain allows for radiating less energy, but the 
circuit power per AP increases with the number of antennas— 
the multiplexing gain may, on the other hand, reduce the circuit 
power spent per UE. Designing cellular networks for high EE 
is thus a non-trivial task, where the throughput, AP density 
and hardware characteristics need to be considered. 

The spatial distribution of UEs is highly heterogeneous in 
practice, which calls for an increasingly heterogeneous and 
complex AP deployment. A promising way to model and ana¬ 
lyze such networks is by using stochastic geometry | |T3 , where 
the AP locations form a realization of a spatial point process, 
typically a Poisson point process (PPP). This approach can 
provide tractable expressions for key performance metrics such 
as the coverage probability and average spectral efficiency 
(SE) in the network. A few prior works have also derived EE- 
related performance expressions and showed how these depend 
on the AP and UE densities; for example, HD considered the 
deployment of two types of single-antenna APs, while |14| 
studied the EE when multi-antenna APs serve one UE each. 

In this paper, we consider networks with M-antenna APs, 
each serving K UEs. In contrast to our previous work |^, 
which considered spatially symmetric AP deployments, we use 
stochastic geometry with a given AP density A to model the 
spatial randomness in network deployment. We obtain a lower 
bound on the achievable EE and maximize it analytically with 
respect to M, K, A, and the transmission power p. The result¬ 
ing expressions reveal the fundamental interplay between these 
four design parameters, which are also illustrated numerically. 


II. System Model 

This paper considers a heterogeneous network with outdoor 
APs that provide wireless broadband access to stationary UEs 
in a two-dimensional coverage area. The APs are distributed 
in according to a homogeneous PPP ikA of intensity 
A [APs per m^]. Each AP is equipped with M antennas 
and communicates with K single-antenna UEs, which are 
uniformly distributed in the Voronoi cell that the AP has as 
coverage area; see Pig. [T] The stationarity of the UEs implies 
that the coherence interval is sufficiently large so that the 
overhead required for channel estimation can be reasonably 
neglected. The channel estimation errors are also neglected 
as the distortion noise caused by hardware impairments is 
assumed to be the dominant one G3- These assumptions are 
well-justified in a wireless broadband access scenario where 




Fig. 1. Illustration of the AP positions modeled by the homogeneous PPP 
Each AP has M = 4 antennas and serves K = 3 star-marked UEs. 


each UE is static (e.g., located in a home) and demands 
continuously high data traffic (e.g., for video streaming). 

We consider the downlink and utilize the translation- 
invariance of PPPs to concentrate on a typical UE, which is 
representative for any UE in the network m- It is connected 
to APq G 'I'a and has an arbitrary UE index k. The received 
baseband signal G C at the typical UE is modeled as 

yk = Vl - j Vph^fcWoSo -f ^ j 

\ je'i'A\{APo} / 

+ ek+nk (1) 

where ^ ^ CJV{0,uj~^d~^lM) is the Rayleigh flat-fading 
channel from AP^ G T'a to the typical UE k, di^k is the 
distance between them, a > 2 is the pathloss exponent, and 
w models fixed propagation losses (e.g., wall penetration). We 
let p denote the average RE transmission energy per symbol. 
Si ~ CN{Q, Ik) contains the normalized information symbols 
sent by AP^, and define G as the corresponding 

beamforming matrix with normalized columns. The additive 
receiver noise is modeled by Uk CN(Q, cr^) with being 
the noise energy per symbol, while e and model the 
distortions from hardware impairments and are described next. 

A. Distortion Noise from Hardware Impairments 

Contrary to most other works on energy efficiency in 
wireless networks, the system model in Q includes distortion 
noise from hardware impairments. These distortions are due 
the unavoidable clock drifts in local oscillators, finite-precision 
analog-to-digital converters, non-linearities, finite-order analog 
filters, etc. Such impairments cannot be fully removed HD’ 
but are negligible at low signal-to-noise ratio (SNR). However, 
the impact at high SNR cannot be neglected since it fully 
characterizes the maximal achievable throughput pS] . 

Similar to JT3-i[Tg, we model the transceiver hardware 
impairments as reduction of the received signals by a factor 
s/1 — and replacing it with Gaussian distortion noise that 
carries the removed power: 

efe~CAA(0,eVE F (2) 

V iG-Pa / 

where || • || is the Euclidean norm. The parameter 0 < e < 1 is 
tightly connected to the error vector magnitude (EVM), which 


is a common measure of the transceiver quality. The EVM 
is defined as the ratio between average distortion magnitude 
and the signal magnitude, which becomes « g with 

our notation. Typical EVM values are in the range e < 0.17, 
where smaller EVMs allow for higher SEs |[20|. 


B. Average Achievable Spectral Efficiency 


The SE of a communication link is strictly upper bounded 
by the channel capacity. In contrast, there is no strict upper 
bound on the EE metric of a communication link; showed 
that one can achieve an unbounded EE in massive MIMO if the 
circuit power is neglected—and only time will tell how small 
the circuit power becomes in future hardware. In other words, 
one can only define achievable lower bounds on the EE of a 
network. In this paper we derive a lower bound that is tractable 
for analytic optimization, and for that we need a closed-form 
SE expression. We consider zero-forcing beamforming (ZEBE) 
since it gives the highest EE in the single-cell analysis of 
The typical ergodic achievable SE with ZEBE i^ 
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where the expectation is with respect to channel fading for a 
given T'a. We now take the average of (|^ over the PER This 
is hard to do in closed form, but we can obtain a lower bound. 


Proposition 1. An achievable lower bound on the average SE 
[bit/symbol/user] of the network with ZEBE (M > K + 1) is 


(4) 


Proof: The proof is based on taking the average of ([^ 
with respect to the channel realizations and the PPP T'a. and 
then computing an achievable lower bound by using Jensen’s 
inequality as E{log 2 (l -t- 4)} > log 2 (l -f g^). 

Computing the expectation with respect to the chan¬ 


nel realizations yields E{ 
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due 

and 


^ZF|2 

to the properties of complex Wishart matrices 
E{l|hffcWPP \d,,k} = Kuj-^d-f for i ^ 0 since thTzPBP 
at APq is independent of channel realizations in other cells. 
The distance to the serving BS is do.fc Rayleigh(- 


zttA ^ 


|23|. We keep fixed and compute the expectation with 


respect to the interfering APs (which are further away) as 
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123 Proposition 2.13]. Einally, we make use of EjdQ;.} = 
for V > —2 (e.g., v = 2 and z/ = a) to get Q. ■ 

The lower bound on the average SE expression in Propo¬ 
sition [T] is used in the remainder of this paper to define an 
achievable EE and to optimize this metric. The tightness of 
the lower bound in Proposition [T] is evaluated in Section |V] 


'The SE in is obtained by canceling intra-cell interference using ZFBF 
(see |21| Sec. 3.4] for details), treating inter-cell interference as worst-case 
Gaussian noise in the decoding, and using the distortion noise variance in 




















III. Metric & Problem Formulation 


EE(p, A, M, K) 


The key performance metric in this paper is the achiev¬ 
able downlink energy efficiency (EE) [bit/Joule]. It is de¬ 
fined as the ratio between the area spectral efficiency (ASE) 
[bit/symbol/m^] and the area energy consumption (AEC) 
[Joule/symbol/m^]. Eor a given SE the ASE is defined as 

ase = aa:se. (5) 


The overall AEC accounts for radiated signal energy, dissipa¬ 
tion in circuits, digital signal processing, backhaul signaling, 
and overhead such as cooling. These are all non-negligible 
parts of the energy consumption of practical systems We 
utilize the generic power consumption model from @ and let 
the AEC take the following general form: 

AEC = A +Co+ CiK -f VoM + V^MK^ + ASE 

( 6 ) 

where rj G (0,1] is the efficiency of the RE amplifiers at the 
AP and we recall that p is the average RE transmission energy 
per symbol per active UE. The term Cq is the static energy 
consumption of an AP, while CiK and VqM are the terms 
that scale linearly with the number of active UEs and number 
of antennas at the AP. The higher-order term ViMK accounts 
for the cost of signal processing, in particular, computing the 
product 'WiSi for each information vector and other matrix 
operations]^ The energy consumed by the coding and decoding 
process (and also backhaul signaling) is proportional to the 
ASE with proportionality coefficient A. More details on this 
energy consumption model can be found in ||^ Section IV]. 
In the analysis, the parameters A, rj, Cj, and Vj are fixed at 
arbitrary positive values. Example values are given in Table 
In summary, the EE is defined as using (j^ and (j^. 


A. Problem Formulation 


We use the lower bound in Proposition to facilitate ana¬ 
lytical optimization of the EE, while numerical optimization is 
done in Section |V] Optimization of a lower bound guarantees 
achievable results. The considered lower bound on the EE is 


EE=- 


AiTSE 


+C 0 +C 1 K + VoM + ViMK + AKSFA 


■ (7) 


Note that the AP density A appears as a multiplicative factor 
in both the numerator and the denominator of 0 and thus 
cancels out. However, the EE still depends on A as it also 
appears in the lower bound SE in Q. 

_As seen from 0 , we can regard the EE as a function 

EE(/ 9 , a, M, K) of the downlink transmission power per UE 
ip), the AP density (A), the number of antennas per AP 
(M), and the number of active UEs per AP (AT). All other 
system parameters are assumed to be fixed in the analysis. 
To guarantee reasonable user performance, the following EE 
optimization problem is considered: 


^The precoding matrix needs to be computed once per coherence 
interval with a complexity proportional to MK^. Since the coherence interval 
is assumed to be very long herein, we neglect this part in the paper. See 
for a full framework that takes also the precoding computation into account. 


maximize 

p>0, 0<A<A„ax 

M,R’GZ+ 

subject to SE = 7 , M > K + I, 


( 8 ) 


where 7 is an SE level that is guaranteed on average to UEs in 
the network]^ Note that p can be any positive number, while 
A is a positive number upper bounded by Amaxi this is the 
highest average AP density that can be physically deployed. 
M and K belong to the set Z+ of strictly positive integers. The 
EE optimization problem in ([^ is solved in the next section. 


IV. Optimizing the Energy Efeiciency 


In this section, we solve the EE optimization problem in ([^. 
While doing so, we will also derive expressions that reveal the 
fundamental interplay between the design parameters. 


A. Optimizing the AP Density 

We begin the optimization by considering the AP density A 
(in AP/m^), when the other parameters are fixed. 

Theorem 1. Consider problem 0 for given values on M, K, 
and p. If the problem is feasible, the EE metric is monotoni- 
cally increasing in A and thus maximized at A* = Amax- 


Proof: The AP density A only appears in 0 as part of 
the SE expression in 0. The EE is an increasing function of 
the SE, which is in turn monotonically increasing in A. Hence, 
the EE is maximized at its highest value Amax- ■ 

This theorem validates the intuition that from an EE per¬ 
spective it is preferable to have as high AP density as possible 
(since the signal and interference power increase accordingly). 
We should keep in mind the assumption that every AP has at 
least one UE. Hence, we want to make the cells as small as 
possible while keeping all APs active. Even under future very 
high user densities the inter-AP distances will be at the order 
of meters, thus making the maximal average AP density Amax 
a finite number (e.g., 0.1 AP/m^). Moreover, this paper uses a 
non-line-of-sight channel model, which might not be valid at 
very short distances. Hence, we treat A as a fixed parameter in 
the range 0 < A < Amax < cc in the remainder of the paper. 


Remark 1. Observe that even if we tef A —> 00 (and 
consequently let p —>■ 0), the EE has the finite upper limit 

Ky 

- - - (9) 

(Co -f CiK + VqM + ViMK + AK-/^ 


because the transmission power term goes away in the EE 
expression, while the circuit power consumption remains. 
Hence, smaller cells will only bring EE improvements till the 
point where the transmission power becomes negligible and 
then higher cell density only brings marginal improvements. 


B. Optimizing the Transmission Power 

Next, we find the optimal transmission power per UE: p*. 
The following theorem shows that the SE constraint in 0 can 
be eliminated by selecting p appropriately. 


^By removing interference and noise, it can be shown that is only 
feasible for 0 < 7 < —2 log 2 (e), thus we assume that 7 lies in this interval. 




Theorem 2. For any values on X, M, and K, the SE constraint 


in (Hj is satisfied by 
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( 10 ) 


if problem is feasible. The problem is infeasible whenever 
p* is negative (i.e., when the denominator of is negative). 


Proof: Follows by solving the SE constraint for p. ■ 
This theorem gives the relationship between p* and other 
system parameters. The optimal transmission power is in¬ 
versely proportional to the AP density as due to shorter 

pathlosses when A increases. It is a decreasing function of the 
number of antennas, M, due to the array gain from coherent 
beamforming; the relationship is as M~^ when M is large. 
Finally, p* increases with K since K makes the denominator 
of ( [Tol l smaller; this is explained by the will to operate at 
higher SNRs when the inter-cell interference grows stronger. 

By substituting p* from Theorem [^ into i) and taking A 
as a constant, the EE optimization problem is reduced to 

. . Kj 

maximize y ^ 
M.ifez+ 

(2'y — l) 2K 

subject to {M-K)> _ 27 ^ 2 ) 

C. Optimizing the Number of AP Antennas 

Next, we find the optimal number of AP antennas, M, when 
p* is given and the other system parameters are fixed. 


Theorem 3. For any given values on A and K, the EE metric 
in is maximized by 

2K{2^-1) I 2'r-l 

(a - 2) (1 - 2T'£2) + Y (^A) t (Po+i?i iT) ■ 

( 12 ) 

If M* is not an integer, then the optimum is attained at either 
the closest smaller or larger integer. 

Proof: Maximizing the EE metric in GD with respect to 
a real-valued M is identical to minimizing the denominator. 
This is a convex problem and ( [T^ is the point where the 
first derivative is zero. The convexity implies that the optimal 
integer M is one of the two closest integers to M* ||^. ■ 

This theorem shows how the optimal number of antennas 
per AP depends on the other system parameters. In particular, 
we see that M* increases roughly linearly with the number 
of UEs. In contrast, fewer antennas should be used if the AP 
density increases (i.e., smaller cells) or when the circuit power 
parameters Vq and Vi grow so that it becomes more costly 
to have additional antennas turned on. Finally, more antennas 
are needed when the SE constraint 7 is increased. 


Theorem 4. For any given values on A and = ^ > 1, the 
EE metric in o is maximized by 


K* 
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If K* is not an integer, then the optimum is attained at either 
the closest smaller or larger integer. 


Proof: Maximizing the EE metric in GD with respect to 
a real-valued K is identical to minimizing EE , which is a 
convex problem. Finding the AT > 0 where the first derivative 
is zero gives a quadratic polynomial equation, from which ( [T3| ) 
is the positive root. The convexity implies that the optimal 
integer-valued K is one of the two closest integers. ■ 

This theorem shows how the number of active UEs depends 
on the other system parameters. We notice that K* is an 
increasing function of the static energy consumption Cq and a 
decreasing function of the circuit coefficient Vi that represents 
the power consumed by signal processing (which behaves as 
T>iMK). Interestingly, Ci and T>q have no impact on K*. 
Higher cell density (i.e., small cells) implies fewer UEs. 


E. Alternating Optimization 

To summarize, we first noted in Theorem [T] that the AP 
density A should be as large as it physically and practically 
can. Next, we eliminated the transmission power p from by 
finding its optimal value in Theorem [^ Then, Theorems and 
[^showed how to optimize the EE separately with respect to M 
and K. However, we would like to solve the original problem 
(j^ jointly with respect to all the parameters. We propose the 
following alternating optimization algorithm; 

1) Fix A at its largest possible value (e.g., Amax)- 

2) Assume that an initial feasible point {p, M, K) is given; 

3) Optimize K by using Theorem (and update p*); 

4) Optimize M by using Theorem (and update p*); 

5) Repeat 3) - 4) until convergence is achieved. 

This algorithm converges since the EE has a finite upper 
bound and is non-decreasing in each iteration. The mixed- 
integer nature of makes it hard to further assess the 
convergence, but we have the following result. 


Lemma 1. Suppose ([^ is relaxed by treating M and K as 
real-valued. The alternating optimization algorithm converges 
to the global optimum of this relaxed EE optimization problem. 


Proof: By treating K and /3 = ^ as the real-valued 

optimization variables, it is easy to show that EE in ( [TT] i is a 
jointly convex optimization problem. It then follows from |24 


Prop. 6 ] that the alternating optimization algorithm converges 
to the global optimum of the relaxed EE problem. ■ 


D. Optimizing the Number of UEs 

Finally, we find the optimal number of active UEs per AP. 
Since we have the constraint M > K 1 in (j^, a tractable 
optimization of K requires that also M is changed. Hence, 
we let /? = ^ be fixed and optimize M and K jointly. 


V. Simulation Results 

In this section, we provide numerical results that illustrate 


the behaviors that were derived analytically in Section IV A 


number of hardware and propagation parameters appear in the 
ASE and AEC models, and these need to be selected in the 


















TABLE I 

Simulation Parameters 


Parameter 

Symbol 

Value 

Pathloss exponent 

OL 

3.76 

Fixed propagation loss 

UJ 

35 dB 

Power amplifier efficiency 

V 

0.39 

Level of hardware impairments 

e 

0.05 

Symbol time 

s 

2 W [s/symbol] 

Coding/decoding/backhaul 

A 

1.15[J/Gbit] 

Static energy consumption 

Co 

lOW-S [J/symbol] 

Circuit energy per active UE 

Cl 

O.IW-S [J/symbol] 

Circuit energy per AP antenna 

Do 

1 W ■ S [J/symbol] 

Signal processing coefficient 

Vi 

1.56 ■ 10“^'^ [J/symbol] 

Noise variance 


10~^^ [J/symbol] 
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Smaller cells = 
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Smaller cells = 

Large EE improvement 

-) = 4 bit/symbol (upper) 
- 0 - -7 = 4 bit/symbol (lower) 
—#—7 = 5 bit/symbol (upper) 
—©—= 5 bit/symbol (lower) 
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Fig. 2. Energy efficiency (in Mbit/Joule) as a function of the AP density in 
AP/km^, for different SE constraints. The number of antennas, number of 
UEs, and transmission power are optimized to yield maximum EE. 


simulations. The parameter values are summarized in Table |I] 
and builds upon a similar list of parameter values in ||^. 
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Fig. 3. Energy efficiency (in Mbit/Joule) for 7 = 3 and A= 10“^ AP/m^ = 
10^ AP/km^. The global optimum is star-marked, while the convergence of 
the alternating algorithm from Section [IV-E|is indicated with circles. 


parameters. Fig. shows the EE lower bound as a function 
of the number of AP antennas (M) and number of UEs 
(K), for optimized transmission powers. The strong inter-cell 
interference creates a need to amplify the desired signals by an 
array gain from coherent beamforming; many more antennas 
than UEs are required to support the prescribed SE. The global 
EE maximum is achieved by {M*,K*) = (193,21), and a 
total transmission power of ^ = 424 mW. This type of 

configurations is known as massive MIMO S-O- 
The alternating optimization algorithm from Section [TV-E| is 
also considered in Eig.|^ It was initiated at (M, K) = (10,1) 
and converged in three iterations to {M,K) = (183,20) with 
an EE of 5.71 Mbit/Joule. The 0.2% deviation from the global 
optimum is due to rounding effects, since optimal converge is 
only guaranteed for real-valued M and K (see Lemma [T]|. 


A. Optimizing Energy Efficiency 

There are four variables in the EE optimization problem of 
but we noticed that the AP density A should be as large as 
it can practically be. We therefore study the EE as a function 
of A. Fig. 1^ shows the EE for an optimal choice of M, K, and 
p, for different average SE constraints: 7 G {4, 5} bit/symbol. 
Both the lower bound on the SE in Proposition [T] and an upper 
bound based on Monte-Carlo simulations are shown. There is 
a gap between these curves, but the behavior is the same—it 
is basically only a scaling difference. As expected, increasing 
the AP density will increase the EE monotonically since the 
propagation losses reduce. However, when the AP density goes 
beyond A = 10“‘*AP/m^ = lO^AP/km^ the improvements 
are really minor. This is because the transmission power is now 
negligible as compared to the circuit power. An AP density 
of A = 10“‘*AP/m^ corresponds, roughly speaking, to an 
average inter-AP distance of 100 meters, which is small but 
not remarkably small. The behavior is the same for both SE 
constraints, but 4 bit/symbol gives the highest EE values. 

Next, we fix the AP density to A = 10“^AP/m^ and 
the SE to 3 bit/symbol, while studying the impact of other 


B. Optimization under Fixed UE Density 

Suppose the APs are deployed to match a certain density 
of UEs: p. UEs per m^. To accommodate all these UEs, the 
network should be designed under the additional constraint 

p = KX. (14) 

We now study how this UE density determines the AP density. 
Euture densities from 10^ UEs per km^ (in rural areas) to 10® 
UEs per km^ (in shopping malls) have been predicted in the 
METIS project | p5| , and are used as reference points herein. 

Fig. 1^ shows the EE has a function of the UE density for 
7 = 3 bit/symbol, while Eig. shows the corresponding AP 
density. The design parameters M, K, A, and p are optimized 
as in (|^ but with the additional constraint in ( [l4| l. Two 
reference cases are also shown: Single-user transmission with 
(M, K) = (10,1); and massive multi-user MIMO transmission 
with (M, K) = (195, 20). Only the AP density and transmis¬ 
sion power were optimized for EE in the two reference cases. 

Several important observations can be made. Firstly, almost 
the same EE can be maintained irrespective of the UE density. 
This is mainly achieved by scaling the AP density linearly 




























Fig. 4. Energy efficiency (in Mbit/Joule) as a function of the UE density 
fi = KX. The EE is optimized with respect to {M, K, X, p), or only with 
respect to (A, p) for given M and K. 


with the UE density, while basically the same number of AP 
antennas and UEs per AP are used. Secondly, the hxed massive 
MIMO conhguration {M,K) = (195,20) achieves nearly the 
optimized EE, which is why we chose this reference case. In 
contrast, single-user transmission in each cell leads to much 
lower EE. We stress that the optimal conhguration (and the 
hxed massive MIMO case) serves the UEs with the same SE 
of 3 bit/symbol but using 20 x lower AP density—this is likely 
to be a key property for cost-efficient AP deployments. 

VI. Conclusions 

We have designed cellular networks for energy efficiency. 
This was formulated as an optimization problem by using 
stochastic geometry, a lower bound on the SE, and a state- 
of-the-art power consumption model. The variables were AP 
density, number of antennas and UEs per AP, and the transmis¬ 
sion power. The results show that the EE increases with the AP 
density, but the positive effect saturates when the circuit power 
dominates over the transmission power. A further leap in EE is 
achieved by enabling massive MIMO transmission, by having 
many AP antennas and spatially multiplex many UEs/cell. The 
gains come from intra-cell interference suppression and by 
sharing the circuit power cost between these UEs. The analysis 
focused on providing wireless broadband access to static users, 
while future work will consider fast-fading channels where the 
channel estimation overhead and errors cannot be neglected. 
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